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1. Introduction 

Searches for rare-decay processes have played an important role in the development of the 
Standard Model (SM). The absence of flavor-changing neutral current (FCNC) in K decays im- 
plied the existence of charm quark, and the observations of — and Bg — Bs mixings (FCNC 
processes) signaled the very large top-quark mass. The study of the FCNC has been focused on 
rare decays involving transitions such as 5 — > dl^l^ , s ^ dvv, b ^ sy, and b ^ l^l^. 

In contrast with the K and B FCNC processes, the D meson FCNC transitions are mediated 
by light down-quark sector, which implies an efficient Glashow-Iliopoulos-Maiani (IM) cancella- 
tion. In this case, D" mixing, as well as FCNC decays are expected to be very small in the SM 
short-distance contributions. Many extensions of the SM may enhance the mixing rate and FCNC 
processes in the short-distance of the D system by orders of magnitudes. Unfortunately, the new 
physics contributions are diluted by large long-distance contributions which are likely to dominate 
over the SM short-distance effects. For example, the charm radiative decays D — > /jy are com- 
pletely dominated by hadronic uncertainty, where h denotes a light hadron final state |[T]]. However, 
semileponic decays such as c ^ ul^l^ may be used to constrain the new physics since one can look 
at the new physics contributions to the whole kinematics to which the new physics may contribute 
the region away from the resonance-dominated (long-distance dominated) region Although 
the long-distance effects dominate the mixing rate, the recent experimental bounds can still be 
used to constrain the new physics space. Purely leptonic flavor-violating (LFV) modes such as 
—>■ pL^e^ and D —>■ hpL^e^ , as well as lepton-number violated (LNV) modes D+ h^e^e^ , 
are completely not allowed in the SM, they are "smoking gun" for new physics searches. In this 
paper, we report the recent measurements of rare and forbidden charm decays from CLEO-c, and 
discuss the prospect of these decays at the BES-III experiment. A comparison between B factories 
and T-charm factories are also presented. 

2. CLEO-c and BES-III for charm flavor physics 

The Cornell Electron Storage Ring (CESR) had been upgraded to CESR-c with the installa- 
tion of 12 wiggler magnets to increase damping at low energy. The CLEO-c detector is minimal 
modification of the well understood CLEO-III detector. It is the first modern detector to operate at 
charm threshold. The CLEO-c has accumulated a total of 281 pb^^ at i/a(3770) peak (1.8 million 
DD pairs) and 200 pb^^ at ^/s = 4170 MeV for physics. CLEO-c expects to take data until 
April 2008 and will approximately triple each data set by that time [Q]. 

The BES-III detector is designed for the e+e^ collider running at the T-Charm energy region, 
called BEPCII, which is currently under construction at IHEP, Beijing, PR. China The accel- 
erator has two storage rings with a circumstance of 224 m, one for electron and one for positron, 
each with 93 bunches spaced by 8 ns. The total current of the beam is 0.93 amp, and the crossing 
angle of two beams is designed to be 22 mrad. The peak luminosity is expected to be lO^^cm^^s^' 
at the beam energy of 1.89 GeV, the bunch length is estimated to be 1.5 cm and the energy spread 
will be 5.16 x I0~^. At this moment, the LINAC has been installed and successfully tested, all 
the specifications are satisfied. The storage rings have installed, and will be commissioned for 
synchrotron radiation run by the end of the year. 



2 



Rare and forbidden decay of charm mesons at CLEO-c and BES-III 



Hai-Bo Li 



The BES-III detector consists of a He-based small cell drift chamber, Time-Of-Fhght (TOP) 
counters for PID, a CsI(TI) crystal calorimeter, a solenoid super-conducting magnet with a field 
of 1 Tesla and the magnet yoke interleaved with Resistive Plate Chambers (RPC) counters as the 
muon chamber. The construction is expected to be completed in the middle of 2007. Photon energy 
resolution isAE/E = 2.5% atEy =1.0 GeV. The momentum resolution is Op/p = 0.5% at = 1.0 
GeV/c, and the dE/dx resolution for hadron tracks is about 6%. The time resolution of TOP is 
about 100 ps, combining energy loss (dE/dx) measurement in the draft chamber, give 10 a KJn 
resolution across the typical kinematic range. 



Table 1: r-Charm productions at BEPC-II in one year's running(lO^s). 





Central-of-Mass energy 


Luminosity 


#Events 


Data Sample 


(MeV) 


(lO^^cm^^s-i) 


per year 


J/v 


3097 


0.6 


10 X 10*^ 




3670 


1.0 


12 X 10^ 


Y{2S) 


3686 


1.0 


3.0 X 10*^ 


dOd° 


3770 


1.0 


18 X 10^ 


D+D 


3770 


1.0 


14 X 10*^ 




4030 


0.6 


1.0 X 10^ 




4170 


0.6 


2.0 X 10^ 



The BES-III can accumulate 10 x lO"^ J/\if, 3 x lO'^ v^(2S) , 30 million DD or 2 million D^D^- 
pairs per running year as listed in Table [D, respectively, when it is turned to run at resonances in 
2008 ||^. Coupled with what is available at CLEO-c, the BES-III will make it possible for the first 
time to study in detail the light hadron spectroscopy in the decays of charmonium states and the 
charmed mesons. In addition, about 30 million DD pairs will be collected at BES-III in one year at 
yif{311Q) peak. With modern techniques and unprecedented high statistical data sample. Searching 
for rare D, Charmonium and tau decays will be possible, such as LEV and LNV decays, invisible 
decays. 

3. Rare charm decays from CLEO-c and prospect at BES-III 

With 1.6 million charged D at CLEO-c, they have searched for the FCNC decays D+ — > 
TT+e+e^ and D+ e+e^, and the LNV decays D+ n^e^e^ and D+ —>■ K^e^e^ (charge- 
conjugate modes are implicit throughout this paper) [^]. For each candidate decay of the form 
D+ —>■ h'^e^e^, where h is either n or K, the energy difference AE = E^and — Eheam is computed and 
the beam-constrained mass difference Ambc = \J E^eam ~ \Pcand\'^ where Ecand and Pcand are 

the measured energy and momentum of the h'^e^e^ candidate, Ebeam is the beam energy, and m[)+ 
is the niminal mass of the D+ meson. Event with D+ candidates satisfying —30 MeV < Am^c < 30 
MeV and -100 MeV < AE < 100 MeV are selected for further study. Within this region, the signal 
box is defined as -5 MeV < Arrive < 5 MeV and -20 MeV < AE < 20 MeV, corresponding to ±3a 
in each variable, as determined by MC simulation. The remainder of the candidate sample was 
used to assess backgrounds. 
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Figure 1: Scatter plots of Amh^ vs. AE obtained from data for each decay mode. The signal region, defined 
by -20 MeV < AS < 20 MeV and -5 MeV < Am^c < 5 MeV, is shown as a box. The two contours for 
each mode enclose regions determined with signal MC simulations to contain 50% and 85% of signal events, 
respectively. 



Candidates for the decay of the long-distance decay D+ — > 7r+(/) — > n^e^e^ are selected based 
on the mass squared of the final-state e^e^ (equal to the of the decay), with 0.9973 GeV^ 
< nr^+g- < 1.0813 GeV^ defining the -resonant region. This region is used both to veto the long- 
distance D+ ^ 071+ — > TT+e+e^ contribution and to measure its branching fraction. Backgrounds 
in the D+ K^e^e^ candidate sample arise from both DD and non-DD sources. A detailed 
analysis of backgrounds can be found in Ref. [0]. 

CLEO-c results show no evidence for signals, the 90% confidence level upper limits are set as 
below: 

SS{D+ K+e+e-) < 7.4 x 10"^ 
^{D+ K-e+e+) < 3.6 x 10"^ 
^(D+ ^ K+e+e-) < 6.2 x 10"^ 
^{D+ ^K-e+e+) < 4.5 x 10^^ 

these results for these dielectron modes are significantly more restrictive than previous limits, and 
reflect sensitivity comparable to the searches for dimuon modes [16]. Due to the dominance of 
long-distance effects in FCNC modes, they separately measure the branching fraction of the reso- 
nant decay D+ 71+0 K+e+e^, obtaining ^{D+ 071+ Ti+e+e^) = (2.7;j-^±0.2) x 10"^. 
This is consistent with the product of known world average [[T^ branching fractions, ^{D+ — > 
4>n+ n+e+e-) = ^{D+ <pn+) x ^{(j) e+e') = [(6.2 ± 0.6) x lO^^] x [(2.98 ± 0.04) x 
10-4] = (1.9 ±0.2) X 10-^ 

With 20 fb^^ data at i/a(3770) peak with the BES-III detector, the sensitivity of the measure- 
ments of the rare charm decays are summarized for D+ and D", respectively, in Tables ^ and |3[ For 
most of them, the sensitivities can be at the order of 10^^. The current best experimental hmits are 
also listed in the tables ^ and |3| 



4. - D° Mixing, CP Violation at BES-III 

With the design luminosity of lO^^cm^^s^^, BES-III will have the opportunity to probe for 
the possible new physics which may enter up-type-quark decays. It includes searches for charm 
mixing, CP violation and rare charm decays. The BES-III charm physics program also includes a 
variety of measurements that will improve the determination of 03 // from B-factory experiments. 
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Table 2: Current and projected 90%-CL upper limits on rare D+ decay modes at BES-111 with 20 fb data 
at V/(3770) peak. We assume the selection efficiencies for all modes are 35%. 





Reference 


Best Upper 


BES-III 


Mode 


Experiment 


limits(10"^) 


(xlO-^) 




CLEO-c 


7.4 


5.6 




FOCUS [|] 


8.8 


8.7 




E791 [0] 


34 


8.3 




CLEO-c [0] 


3.6 


5.6 




FOCUS [|] 


4.8 


8.7 




E791 [|] 


50 


5.9 


K+e+e- 


CLEO-c [0] 


6.2 


6.7 




FOCUS [|] 


9.2 


10.5 




E791 [|] 


68 


8.3 


K-e+e+ 


CLEO-c 


4.5 


6.7 


K 


FOCUS [|] 


13 


10.4 


K-n+e+ 


E687 [|10|| 


130 


8.3 



Table 3: Cun-ent and projected 90%-CL upper limits on rare decay modes at BES-111 with 20 fb ^ data 
at 1/(3770) peak. 



Mode 



Reference 
Experiment 



Best Upper 



limitsClO"*^) 



BES-III 

(xlO-^) 



77 



Ksix+ix- 
Kse+e- 



CLEO [112p 
BaBar |[T|] 
BaBar |[T|] 
BaBar 
E653 [HI] 
CLEO [HI 



CLEO [jlj 
E653 [O] 



CLEO [13] 



CLEO [13] 



CLEO [13] 



CLEO [13] 



CLEO [13] 



28 


5.0 


1.3 


17 


0.81 


4.3 


1.2 


2.4 


180 


12.3 


86 


9.7 


45 


7.9 


260 


10.6 


100 


9.6 


110 


7.5 


530 


10 


100 


10 


110 


10 
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The total number of charm mesons accumulated at BES-III will be much smaller than that at the 
B-factories which are about 500 million DD pairs for each of them. However, the quantum correla- 
tions in the i//^(3770) DD system will provide a unique laboratory in which to study charm [jl7|]. 

4.1 D"-D° Mixing 

£)0 _£)0 juixing within the SM are highly suppressed due to GIM mechanism, thus, at BES-III, 
searches for neutral charm mixing and CP violation in the charm decays may be essential to search 
for some intriguing signals due to new physics. 

The time evolution of D" — D^ system, assuming no CP violation in mixing, is governed 
by four parameters: x = Am/F and y = AT /IF which are the mass and width differences of D 
meson mass eigenstates and characterize the mixing matrix, 5 the relative strong phase between 
the Cabibbo favor (CF) and the doubly-Cabibbo suppressed (DCS) amplitudes and Rjj the DCS 
decay rate relative to the CF decay rate. The mixing rate Rm is defined as ^{x^ +y'^). The SM 
based predictions for x and y, as well as a variety of non-SM expectations, span several orders of 
magnitude [|], |18|, 19] which is x ~ 3^ ~ 10^^. Presently, experimental information about charm 



mixing parameters x and y comes from the time-dependent analyses. 

At the B factories, the wrong-sign (WS) process, D —>■ K^n^, is used to extract the D" mixing 
parameters by fitting the time-dependent decay rates. The WS process can proceed either through 
direct doubly-Cabibbo-suppressed (DCS) decay or through mixing followed by the right-sign (RS) 
Cabibbo favored (CF) decay D" ^ D° ^ K^n^ . The two decays can be distinguished by the 
decay-time distribution. For |x|, j| ^ 1, and assuming negligible CPV , the decay-time distribution 
for D** K^n^ can be expressed as 



Rwsit) =e-'''{RD + VREy'{rt) + '^^^{rtf), (4.1) 



where x' = xcos5 +3'sin5 and y' = — xsin5 +ycos5, and 5 is the strong phase between the DCS and 
CF amplitudes. Recently, a time-dependent analysis in D ^ Kn has been performed based on 384 
fb"^ luminosity at T(45) by BaBar experiment [20]. By assuming CP conservation, they obtained 



the following neutral D mixing results 

Rd = (3.03±0.16±0.10) X 10-^ 
x' = (-0.22 ±0.30 ±0.21) X 10"^ 

y = (9.7±4.4±3.1) X 10"^ (4.2) 

The result is inconsistent with the no-mixing hypothesis with a significance of 3.9 standard devia- 
tions. The results from BaBar and Belle are in agreement within 2 standard deviation on the exact 
analysis of y' measurement by using D Kn as listed in Table |^ 

At i/a(3770) peak, to extract the mixing parameter y, one can make use of rates for exclusive 
D"D° combination, where both the D° final states are specified (known as double tags or DT), as 
well as inclusive rates, where either the D" or D" is identified and the other D° decays generically 



(known as single tags or ST) ^^]\. With the DT tag technique p4] , Pq ], one can fully consider 
the quantum correlation in C = — 1 and C = ±1 D^D^ pairs produced in the reaction e^e^ — > 
D^D^{mt^) and e+e- D^Df^y{nn^) |E|,|2^, |T7|], respectively. 
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Table 4: Experimental results used in the paper. Only one error is quoted, we have combined in quadrature 
statistical and systematic contributions. 



Parameter 



y 

Rd 

ycp 

X 

y 



BaBar (xlO^-^) 



-0.22 ±0.37 [|^: 

9.7 ±5.4 1^] 
3.03±0.19 [Eg 



BeIle(xlO-^) 



+0.21 
-0.23 



[221 



0.18 



o.6tti [© 



3.64±0.17 [22| 



13.1±4.1 [21] 



8.0±3.4[g|] 
3.3±2.8 [|23|1 



Technique 



Kn 
Kn 
Kn 

K+K , K+K 



For the ST, in the limit of CP conservation, the rate of D° decays into a CP eigenstate is given 



as [17]: 



r/,^r(D0^/,) = 2A2 [l-r7j]^ 



(4.3) 



where is a CP eigenstate with eigenvalue Tj = ±1, and Ay^^ = \{fri\J^\D^)\ is the real-valued 
decay amplitude. 



For the DT case, Gronau et. al. [ |280 and Xing [g9|] have considered time-integrated decays 
into correlated pairs of states, including the effects of non-zero final state phase difference. As 
discussed in Ref. [Q, the rate of {D^D^f=-'^ {l^X){fj^) is described as [||]: 



r/;/,=r[(/±X)(/,)] =A2,^A}^(l+y 



(4.4) 



where Aj±x = \ {l^X\M'\D^) \ is real-valued amplitude for semileptonic decays, here, we neglect 
term since y <^\. 

For C = — 1 initial D^'D" state, y can be expressed in term of the ratios of DT rates and the 



double ratios of ST rates to DT rates [jlTp: 

1 fnuTf_ 



(4.5) 



For a small y, its error, A(3;), is approximately 1 / \/Ni±x, where Ni±x is the total number of (/*X) 
events tagged with CP-even and CP-odd eigenstates. The number Ni±x of CP tagged events is 
related to the total number of D^D" pairs A^(D°D°) through Ni±x ^ N{D^D°)[^{D° -^l^+X)x 
^{D^ ^ f±) XEtag] ~ l-5x 10^^A'^(D°D''),here we take the branching ratio-times-efficiency factor 
{^{D^ f±) X Stag) for tagging CP eigenstates is about 1.1% (the total branching ratio into CP 
eigenstates is larger than about 5% [[T6|]). We find 



My) 



±26 



VW(D0D0) 



±0.003. 



(4.6) 



If we take the central value of y from the measurement of ycp at Belle experiment [21], thus, at 
BES-III experiment, with 20 fb~^ data at \^{311Qi) peak, the significance of the measurement of y 
could be around 4.3 a deviation from zero [|3C 
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We can also take advantage of the coherence of the D" mesons produced at the i/a(3770) peak 
to extract the strong phase difference 5 between DCS and CF decay amphtudes. Because the CP 



properties of the final states produced in the decay of the va(3770) are anti-correlated [p6|,|27|], one 
state decaying into a final state with definite CP properties immediately identifies or tags the CP 



properties of the other side. As discussed in Ref. [28], the process of one D" decaying to K 71+ , 
while the other D*' decaying to a CP eigenstate fr^ can be described as 

^ A^a}^ (1+ 27] V^cos5), 

(4.7) 

whereA = |(^^7r+| J^|D")| andAy^ = |(/^|J^|D°)| are the real-valued decay amplitudes, and we 
have neglected the terms in Eq. ( |4.7| ). In order to estimate the total sample of events needed to 
perform a useful measurement of d, one defined [28, I]] an asymmetry 



i K7i;f+ - i KnJL 

r Tr ' ^ 



where rK7i:f± is defined in Eq. {4J_), which is the rates for the 1/^^(3770) — > D*'d" configuration 
to decay into flavor eigenstates and a CP-eigenstates f±. Eq. ( p7^ implies a small asymmetry, 
£/ = 2^/Rj}Cos5. For a small asymmetry, a general result is that its error A^/ is approximately 
1 / yJNx-iz+, where Nx-iz+ is the total number of events tagged with CP-even and CP-odd eigen- 
states. Thus one obtained 

The expected number Njc-n+ of CP-tagged events can be connected to the total number of D^D" 
pairs A^(D°D0) through A^jf-;j+ ^N{D^D^)m{D^ ^ K-K+) x ,^{D^ ^ f±) x w4.2 x IQ-'^NiD^D^) 
here, as in Ref [^, we take the branching ratio-times-efficiency factor SSiD^ — > x ttag = 1.1%. 
With the measured Pd = (3.03 ±0.19) x 10"^ and mijf K-71+) = 3.8% [111], one found [p|] 



±444 

A(cos5) « , (4.10) 

At BES-III, about 72 x 10^ D°D^ pairs can be collected with 4 years' running. If considering both 
K^Ti^ and K~^n~ final states, we thus estimate that one may be able to reach an accuracy of about 
0.04 for cos5. 

At BES-III, the measurement of Rm can be performed unambiguously with the following 



reactions [26]: 



(/) e+e~ 


V^(3770) - 






(a) e^e^ 


V^(3770) - 




{K-e+v){K-e+v 


{Hi) 


-^D D*+ - 




-){7liJK^e-v]). 



(4.11) 
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Reaction (/) in Eq. ( [4.11| ) can be normalized to D^D^ —>■ {K n^){K^n ), the following time- 
integrated ratio is obtained by neglecting CP violation: 

N[{K-K+)iK-K+)] 

N[{K-n^){K^n-)] ^ ~1~ = ^'^^ ^^"^^^ 
For the case of semileptonic decay, as (//) in Eq. ( [4.1 1[ ), we have 



Rm, (4.13) 



In the limit of CP conservation, by combing the measurements of x in D" —>■ Ksun and ycp 
from Belle, one can obtain Rm = (1.18 ±0.6) x 10"*^. With 20fb-' data at BES-III, about 12 
events for the precess D"D° (/r^7r^)(/r*7r^) can be produced. One can observe 3.0 events after 
considering the selection efficiency at BES-III, which could be about 25% for the four charged 
particles. The background contamination due to double particle misidentification is about 0.6 event 



with 20fb-' data at BES-III [31 1. 



4.2 CP Violation in D System 

For the direct CP violation, the SM predictions are as large as 0.1% for decays, and 1% 
level for and Ds decays p2|]. At BES-III, one can also look at the CP violation by exploiting 
the quantum coherence at the i/a(3770). Consider the case where both the and the D" decay into 
CP eigenstates, then the decays \\r{3110) or /i/i are forbidden, where /+ (/_) denotes 

a CP+ eigenstate {CP- eigenstate). This is because CP{f!^f!^) = CP{fi^)CP{fi){-\y = -1, 
while, for the / = 1 i//^(3770) state, CP(i//^(3770)) = +1. Thus the observation of a final state 
such as {K^K^){n^7i^) constitutes evidence of the CP violation. For {K^K^){n^7i^) mode, the 
sensitivity at BES-III is about 1% level. Moreover, all pairs of the CP eigenstates, where both 
eigenstates are even or both are odd, can be summed over for the CP violation measurements at 
BES-III. 



5. Summary 

In summary, the rare charm decays are searched for at CLEO-c. No evidence is found either 
for the rare (FCNC) decays or for the forbidden (LNV) decays of charged D mesons to three-body 
final states with dielectrons. The sensitivities of the rare charm decays are presented by assuming 
20 fb^' at i//^(3770) peak with BES-III detector. The measurements at BES-III for these rare decays 
can be at the level of 10^^. The sensitivities for the neutral D mixing and CP violation are also 
studied, we find that sensitivity for the lifetime difference y can be 0.003 with 20 fb^^ data at 
i/a(3770) peak, and the significance will be 4.3 a if the y is at 1% level. 
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